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Abstract 

In our recent work on concentrated suspensions of uniformly porous colloidal spheres with ex- 
cluded volume interactions, a variety of short-time dynamic properties were calculated, except for 
the rotational self-diffusion coefficient. This missing quantity is included in the present paper. 
Using a precise hydrodynamic force multipole simulation method, the rotational self-diffusion co- 
efficient is evaluated for concentrated suspensions of permeable particles. Results are presented 
for particle volume fractions up to 45%, and for a wide range of permeability values. From the 
simulation results and earlier results for the ffi'st-order virial coefficient, we find that the rota- 
tional self-diffusion coefficient of permeable spheres can be scaled to the corresponding coefficient 
of impermeable particles of the same size. We also show that a similar scaling applies to the 
translational self-diffusion coefficient considered earlier. From the scaling relations, accurate ana- 
lytic approximations for the rotational and translational self-diffusion coefficients in concentrated 
systems are obtained, useful to the experimental analysis of permeable-particle diffusion. The 
simulation results for rotational diffusion of permeable particles are used to show that a general- 
ized Stokes-Einstein-Debye relation between rotational self-diffusion coefficient and high-frequency 
viscosity is not satisfied. 
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I. INTRODUCTION 



The rotational and translational self- diffusion of interacting colloidal and macromolecular 
particles suspended in a low-molecular-weight solvent is the subject of ongoing research 
both experimentally and theoretically |l, 2|- Originally, self-diffusion in dilute systems was 
studied, however the center of interest has shifted since to concentrated dispersions where 
solvent-mediated many-particle hydrodynamic interactions (Ills) are of central importance. 
An example, of biological relevance, is self-diffusion of proteins and other macromolecules 

n 

in the crowded environment of a cell p[. 

Two central quantities quantifying the configuration-averaged influence of His on the 
suspensions dynamics are the concentration-dependent short-time rotational and transla- 
tional self-diffusion coefficients Dr and Dt-, respectively. At zero particle concentration, 
these quantities reduce to the single-particle diffusion coefficients Dq and D^. For solvent- 
impermeable colloidal hard spheres with stick hydrodynamic surface boundary conditions, 
the single-particle coefficients are given by 

" 8^ 
Qirrjoa ' 

with ?7o the Newtonian solvent shear viscosity, Boltzmann constant ks, temperature T and 
hydrodynamic particle radius a. The influence of the His at non-zero concentrations gives 
rise to values for Dr and Dt smaller than their respective values Dq and D^ at infinite 
dilution. The short-time coefficients describe self-diffusion on the time scale t <^ a'^/D^, 
but with t large enough that solvent and particle velocity correlations have decayed. On 
the colloidal short-time scale, the concentration dependence of Dr and Dt is determined by 
averaging the His with the equilibrium particle distribution. 

Self-diffusion coefficients in colloidal suspensions have been determined experimentally by 
a variety of techniques. The mean-squared displacement (MSD) of partially solvent-index 
matched suspensions of colloidal spheres j^, 5| has been measured as a function of time using 
dynamic light scattering (DLS), with Dt determined from the initial slope of the MSD. 

For the vast majority of systems where this specialized index-matching technique is not 



applicable, Df may be inferred, to decent accuracy according to theory and simulation 



from a first cumulant analysis of the scattered light electric field autocorrelation function. 



probed at a scattering wavenumber larger than the peak location of the static structure factor 
where the structure factor attains the value one |9-11|. Translational long-time self-diffusion 
coefficients not considered in the present work can be determined using, e^., forced Rayleigh 



scattering 
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13|, fluorescence recovery after photobleaching (FRAP) [2|], and fluorescence 
correlation spectroscopy jl^ . 

Experimental studies of rotational colloidal self-diffusion are based on techniques which 
can distinguish different particle orientations. Methods which have been used for this pur- 
pose are depolarized dynamic light scattering (DDLS) on optically anisotropic particles jl5 |. 



and nuclear magnetic resonance 



16|. More recently developed techniques at 



a larger variety o: 
polarized FRAP 



icable to 

systems include time-resolved phosphorescence anisotropy |17l . fisj ] and 
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20| measurements. The latter methods have been carried out us- 



ing fiuorophore-labeled colloidal particles. Most of the published experimental results deal 
with self-diffusion properties of monodisperse colloidal systems. However, experimental and 
theoretical work has been also performed on rotational diffusion in colloidal mixtures, iii par- 
ticular for binary systems where one component (the tracer) is very dilute ,[i8[y,i2ij. 

In addition, DDLS measurements of the rotational diffusion of tracer spheres in a poly- 
meric solution have been used to infer viscoelastic pro perties from a frequency-dependent 
generalized Stokes-Einstein-Debye (GSED) relation 22|. 

From a theoretical viewpoint, short-time rotational self- diffusion in monodisperse col- 
loidal systems of non-permeable spheres was studied using lattice-Boltzmann (LB) 
Stokesian dynamics {24 1, and accelerated Stokesian dynamics (ASD) simulations ^]. While 
there has been no theoretical work so far on rotational self-diffusion in dense suspensions of 
porous particles, other transport properties of porous particles have been studied, includ- 
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and to first order in concentration the mean 



ing the high-frequency shear viscosity 
sedimentation velocity |29| . 

In our earlier work on the short-time dynamics of concentrated suspensions of uniformly 
porous particles, a broad spectrum of dynamic properties has been calculated, including the 
hydro dynamic function 



coefficient 



30l | , and the high-frequency-limiting shear viscosity rjoo 3ll |32| . These simula- 



tion studies were amended by the derivation of easy-to-use approximate analytic expressions 



of good accuracy, notably a generalized Saito formula for the shear viscosity 



32l | , and a spher- 



ical annulus model approximation for rjoo 32|, and to first order in concentration also for Dt 
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and Dr 33|]. Additionally, precise values for the first-order virial coefficients of Dr and Dt 



corresponding to two-body His have been obtained 



33|. 



In all these studies on permeable particles, the solvent flow inside the spheres is described 



35| . and the particles are assumed 



by the Debye-Biiche-Brinkman (DBB) equation 
to interact directly by excluded volume (i.e., hard-sphere type) forces. Our simplifying 
particle model is specified by two parameters only, namely the particle volume fraction 
(j) = {A7r/3)na^, where n is the number concentration, and the ratio, x, of the particle 
radius, a, to the hydrodynamic penetration depth, inside a permeable sphere. Large 
(low) values of x correspond to weakly (strongly) permeable particles. Typical values for 
X in permeable-particle systems, such as core-shell particles, are in the range of x ~ 30 or 



larger 



36|. While a specific intra-particle structure is ignored in the model, it is generic 



in the sense that a more complex internal hydrodynamic structure can be approximately 
accounted for in terms of a mean permeability. Porous-particle systems of current interest 



include dendrimers 



37 



391], microgel particles |40l-l42| . core-shell particles |36l. I43l-l45| . and 



star-like polymers of lower functionality 46 1. 

The present work complements our earlier analysis of the short-time dynamics in con- 
centrated suspensions of uniformly permeable spheres by giving simulation results and a 
theoretical analysis of the short-time rotational diffusion coefficient not considered so far 
at non-dilute concentrations. On employing the multipole simulation method of a very 
high accuracy encoded in the hydromultipole program package 48], we calculate the 
short-time rotational self-diffusion coefficient, Dr{x, 0), as a function of x and 0. Our results 
cover the full range of porosities, with the volume fraction extending up to 0.45. In combina- 
tion with recently obtained tabulated values for the first-order virial coefficients of Dr{x, 0) 



and Di 
earlier 



x,6) 33|, and precise hydromultipole simulation results for Dt{x,(f)) obtained 



30|, we show that both Dr{x,(f)) and Dt{x,(j)) can be scaled, in the whole range 



of permeabilities and volume fractions, to the self-diffusion coefficients D^^{(f)) = Dr{oo,(j)) 
and D]^{(j)) = Dt{oo,(j)) of non-permeable hard spheres with stick boundary conditions and 
the same size. From these scaling relations, accurate analytic expressions for Dr{x,(f)) and 
Dt{x,(j)) are obtained. We expect these expressions to be useful in the experimental data 
analysis of diffusion measurement on permeable particle systems. The present simulation 
results for Dr{x,(j)), and known results for ?7oo(a^,0), are used to show the violation of a 
GSED relation between Dr{x,(f)) and T]^{x,(f)), amending our earlier study of similar GSE 
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relations in 



3l|. 



The paper is organized as follows: Sec. |TT] provides the theoretical background on short- 
time self-diffusion of permeable particles. Furthermore, it includes our simulation results for 
rotational self- diffusion. The scaling relations allowing to map permeable to non-permeable 
hard-sphere systems are discussed in Sec. Illli In Sec. \IV\ we complete the scaling relations 
by providing simple expressions for the scaling functions for non-permeable hard spheres. 
We also discuss the special case of non-permeable hard spheres in comparison to earlier 
simulations and experimental work. In Sec. IVj we demonstrate the violation of the GSED 
relation. In our conclusions in Sec. Wl\ we explicitly write convenient expressions for Dr{x, (p) 
and Dt{x, 0) which should prove useful in practical applications. 



II. SHORT-TIME ROTATIONAL SELF-DIFFUSION: THEORY AND RESULTS 

Like in our earlier work on the dynamics of permeable particle systems js], IsolsS], we 
employ a model of uniformly permeable spheres of radius a, dispersed in a Newtonian fluid of 



viscosity rjQ. The low- Reynolds number incompressible flow inside and outside 



described, respectively, by the Stokes j49|, |50| and Debye-Bueche-Brinkman 



;he spheres is 



35| equations 



r]o V'^v{r) - r]o x(r) [v{r) - Ui{r)] - Vp(r) = . 



(3) 



Here, v and p are the fluid velocity and pressure, respectively, and is the hydrodynamic 
penetration depth. The characteristic function, x{^)i is equal to one for the field point 
r inside any of the spheres and zero outside. The skeleton of a particle i, centered at 
Tj, moves rigidly with the local velocity Uj(r) = Uj + cUj x (r — Rj), determined by the 
translational and rotational velocities Uj and a;^, respectively. The fluid velocity and stress 
change continuously across a particle surface. 

The short-time rotational self-diffusion coefficient of a quiescent, isotropic system is given 



in frame-invariant notation by 21 



ksT / 1 
3 \N 



N 



(4) 



where X = {ri, ■ ■ ■ ,rjv} is the configuration of 1 sphere centers, and Tr denotes 

the trace operation. The hydrodynamic mobility tensor, yii^f (X), linearly relates the torque 
acting on a particle i to its rotational velocity, for zero forces and torques exerted on the 



other particles. For the present model system, the average (■ ■ ■ ) is taken over an equilibrium 
distribution of non-overlapping spheres, consistent with the periodic boundary conditions 
used in our simulations. Our numerical calculation of Dr{x, 0) makes use of Eq. (jl]). 

The coefficient Df is a function both of x and 0. At infinite dilution, Eq. reduces to 



51 



52| 



knT 



3 coth X 



X 



(5) 



Note here that DUx) > Dn^^ unless 



X 



oo. 



We have calculated Drix, (t>) to high precision using a hydrodynamic multipole method 



corrected for lubrication 



471 . |48|, |53|, Isj], and encoded in the hydromultipole program 



package extended to permeable spheres. The hydrodynamic particle structure enters into 
the HYDROMULTIPOLE method only through a single-particle friction operator, whose form 
is known for a variety of particle models 5ll453l|. The details of the simulation method 
are given elsewhere The values for Dr presented in the following have been determined 
from equilibrium configuration averages over typically = 256 particles in a periodically 
replicated cubic simulation box, using 100 initial random configurations for each set of 
parameters. This gives a statistical relative error of less than 0.001. In our multipole 
expansion method used for the rotational mobility tensor in Eq. (jll), the multipole order, L, 
was truncated usually at L = 3. To gain high-precision data, extrapolations to L — )■ 8 were 
made, leading to an accuracy in Dr better than 1%. The calculated values for Dr{N) using 
the periodic simulation box with particles are not critically dependent on the system 
size, since Dr{N = oo) — Dr{N) scales with the particle number like 1/A^. This system 
size dependence is similar to that of the high-frequency viscosity, ?7oo(a;,0), of permeable 
particles. The latter was calculated in earlier work 251. l32|. 

Table [T] lists our high-precision simulation results for Dr{x, 0), for volume fractions up to 
(f) = 0.45. Values of the inverse (reduced) penetration depth x are considered from a very 
small value x ~ 5, characteristic of highly permeable particles, up to x = oo characteristic 
of dry particles with stick surface boundary conditions. 
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TABLE I: Simulation results for the normalized short-time rotational self-diffusion coefficient 

Dr{x,(P)/D'Q{x). 



0\x 5 10 20 30 50 100 oo 

0.05 0.995 0.987 0.980 0.977 0.973 0.970 0.967 

0.15 0.983 0.958 0.934 0.922 0.911 0.901 0.888 

0.25 0.968 0.925 0.881 0.860 0.839 0.820 0.796 

0.35 0.951 0.886 0.820 0.788 0.757 0.729 0.690 

0.45 0.932 0.842 0.753 0.711 0.669 0.629 0.576 

III. SCALING SELF-DIFFUSION OF PERMEABLE TO NON-PERMEABLE 
PARTICLES 

From analyzing the numerical data for the rotational self-diffusion coefficient in Table [H 
we have found an interesting scaling of permeable to non-permeable spheres of the same 
size. In addition, we found that a similar scaling is valid for the translational self-diffusion 
coefficient. Therefore, results for both quantities will be given in this section. We start from 
a brief comparison of Dr{x, (p) to -Dt(x, (p). 

The simulation results for Dr{x,(j)) from tableware depicted in the left panel of Fig. [T] 
using symbols. For comparison, the right panel of Fig. [1] shows the corresponding simulation 
results for Dt{x, (p) taken from [8']. For permeable particles, the fluid is allowed to penetrate 
so that the strength of the His is decreasing with increasing permeability, i.e., decreasing 
X. This is the reason for the larger values of Dr and Dt at larger permeabilities. Our 
results show that the effect of His on is weaker than on Dt, i.e., for a given x and 0, 
the reduction of the self-diffusion coefficient relative the infinite dilution value is smaller for 
rotational diffusion. 

The numerical results for Dr{x, 0) and Dt{x, 0) plotted in Fig. [1] have significantly differ- 
ent slopes at small volume fractions 0. On the other hand, these slopes are well-reproduced 
by the first-order virial coefficients, Xr{x) and Xt{x), defined by the following relations, 

= 1 + A.(x)0 + O(0^) (6) 

= 1 + A,(a:)0 + O(0^), (7) 



8 



0.05 0.15 0.25 0.35 0.45 0.05 0.15 0.25 0.35 0.45 



FIG. 1: Rotational (left) and translational (right) self-difFusion coefficients, Dr{x,(j)) / Dq[x) and 
Dt{x, (j))/ Dq{x), as functions of (j), for values of x as indicated. Symbols: simulation results. Solid 
lines: interpolated r.li.s. of Eqs. (|13til4p . 



evaluated in Ref. 33| and listed in Table [TTl The single-particle rotational diffusion coeffi- 
cient^ D^x), has been already given in Eq. (Q, and the translational one has the form given 



m 



35|, 



Dlix) 



1 + 



1 



+ 



X cotha; — 1 



(8) 



TABLE II: First-order virial terms, Xr{x) and Xt{x), of the rotational and translational self-diffusion 



coefficients 





5 


10 


20 


30 


50 


100 


oo 




-0.097 


-0.236 


-0.376 


-0.442 


-0.505 


-0.561 


-0.631 


Xt 


-0.569 


-1.060 


-1.416 


-1.550 


-1.661 


-1.746 


-1.832 



Therefore, the idea is to introduce the following scaling functions, 

1 



Ur{x, 0) 
Ut{x,(j)) 



Dl{x) 
Di{x) 



- 1 

- 1 



\r{x) ' 
1 



(9) 
(10) 



For all values of x > 5, the functions Ur{x,(f)) and Ut{x,(j)) do practically not depend on x, 
i.e. they are permeability- independent. Indeed, as shown in Fig. |2l the curves for Ur and Ut 
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as functions of collapse on the corresponding curves for the non-permeable solid spheres, 
i.e. 

Ur{x,(f>) ~ Ur{00,(f>), (11) 
Ut{x,4>) ^ Ur{00,4>) , (12) 

with a relative error less than 3% for x > 10. 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 



FIG. 2: The functions Ur{x,(j)) and ut{x,(j)) are practically independent of x. Symbols: simulation 
results for the indicated values of x. Solid lines: spline fit interpolations of 1*^(00, (p) and ut{oo, x). 



Therefore, the short-time self- diffusion coefficients in suspensions of permeable particles 
are well approximated by the following expressions, 

^ l + A.(a;)M,(oo,0), (13) 

^ 1 + A,(x)«t(oo,0). (14) 

In Fig. [H the solid, continuous lines are not just mere fits to the simulation data, but 
represent the expressions in Eqs. ( |T3lfT^ . i.e. the outcome of the interesting scaling behavior 
of the short-time self-diffusion of permeable particles. The error made in using Eqs. (JT3]- IT^ 
instead of the precise simulation values, is at most 1% for rotational and 3% for translational 
self-diffusion. To complete the analysis, we need to specify the scaling functions for the non- 
permeable solid spheres. This will be done in the next section. 
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IV. SELF-DIFFUSION COEFFICIENTS OF NON-PERMEABLE SPHERES 



We will now use the existing data for non-permeable hard spheres to construct simple 
approximate expressions for the scaling functions Ur{oo,(f)) and Ut{oo,(f)). 

We start with a comparison between our present simulation results for D^^{(j)) = Dr{oo, 0) 
and D]^^{(j)) = Dt(oo,0) for vanishin^^ permeability, and a selection out of a large body of 



15 



16 



18| ) and theoretical (see, e.g. 



QBE 



55 



published experimental (see, e.g. 
58| ) data on impermeable hard spheres. 

Related to this comparison, we note first that Cichocki et al. have derived precise second- 
order virial expansion results 48 1 



r,hs 



D 



hs 



t.hs 



1-0.6310- 0.7260^ + 0(0'' 



1 - 1.83150- 0.2190^ + 0(0^ 



(15) 
(16) 



for the short-time rotational and translational self-diffusion coefficients of neutral hard 
spheres as functions of 0. 

Regarding rotational diffusion. Fig. [3]^a) shows the comparison of our data with Lattice- 
Boltzmann 23| and ASD j3] simulation results, and DDLS experimental data [15] for op- 
tically anisotropic fluorinated polymer particles. The rotational diffusion coefficient as a 
function of has a concave shape, different from that for Dt which is weakly convex. Our 
simulation data for non-permeable particles agree well with the ASD result. The LB data 
at large are somewhat smaller. The key message conveyed by Fig. [3]^a) is that the second- 
order virial result for D^^(0) in Eq. ( |T5l) describes the simulation and experimental data 
remarkably well for all volume fractions up to the freezing transition value 0.49 {3, |23|, in- 
dicating that higher-order virial coefficients are small or mutually cancel out. Therefore, 
for constructing a simple approximation for Mj,(oo,0) from Eq. (Q, it is sufficient to take 
as Dr{oo,x) the 2nd-order virial expansion in Eq. ( IT5|) . In this way, the rotational scaling 
function is approximated by 



+ 1.151, 



(17) 



with an accuracy of 1.5% or better relative to our simulation data. 



In Fig. [3]^b), we compare our simulation data for D^{(j))/D^ ^ with ASD simulation 



and force multipole calculation results |55|], and with DLS experimental data |6|. The figure 
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(b) 



0.1 0.2 0.3 0.4 



0.5 



0.1 0.2 0.3 0.4 



0.5 



FIG. 3: (a) Rotational and (b) translational self-diffusion coefficients of non-permeable hard spheres 
with stick boundary conditions, as functions of (j). Compared in (a) are our HYDROMULTIPOLE data 



(labeled HYDRO) for D^^"" / Df with SD 7| and LB 



23| simulation results, and DDLS experimental 



data [13] for optically anisotropic particles. In (b), we compare the hydromultipole data for 



551], and DLS 



D]^^ {(j)) / Dq^^ with ASD simulation results [tI, force multipole calculations by Ladd 
experimental data by Segre et al. [6|]. Solid lines: 2nd-order virial expansion results, in (a) according 
to Eq. (fTSj) . and in (b) according to Eq. (fTHl) . 



shows that the translational second-order virial expression in Eq. fll6l) for Df^ noticeably 
underestimates the simulation and experimental data when is larger than 0.3. 

For this reason, we need a more precise expression for Dt{oo, (p) than the 2nd-order virial 
expansion in Eq. (fT6l) . We have found that our simulation data are approximated with 
a 0.4% accuracy by the following expression for the scaling function ut{oo,(f)), defined in 
Eq. 

Ut{oo,(p) ^ + 0.1202 -O.650^ ^^g^ 

The term + 0.120^ follows from the virial expansion in Eq. ( fT6l) . and the coefficient 
of the third order term, —0.65 0^, has been obtained by fitting to the numerical data for 
Ut{oo, 0) - - 0.1202, ^Yie range < < 0.45. 
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V. GENERALIZED STOKES-EINSTEIN-DEBYE RELATION 



We proceed with the discussion of a generahzed short-time GSED relation. Having ob- 
tained in this paper precise numerical data for Dr{x, 0), and taking values of rioo{x, (p) tab- 



ulated in 
relation 



32|, we are in the position to test the validity of the following short-time GSED 

Dr{x,(j)) rioo{x,(j)) ? . 



(19) 



between the rotational self-diffusion coefficient, and the high-frequency viscosity of perme- 
able particles. The validity of generalized Stokes-Einstein relations such as the present one is 
an important issue in microrheological studies where one tries to infer rheological properties 
more easily from diffusion measurements. The GSED relation in Eq. (fT9l) was shown before 
to be violated for suspensions of non-permeable neutral and charged particles 18|, |59(|. Here, 
we ask the same validity question for permeable particle systems. 

In Fig. m the GSED relation is examined for different values of x. If valid, all curves 
should collapse on a single horizontal line of unit height. One notices from the figure that 
the GSED relation is significantly violated for x > 30, and volume fractions > 0.15 
where the particles are significantly correlated. Thus, a rotating particle experiences its 
neighborhood not just as a structureless medium characterized by the viscosity rioo{x,4>). 




0.05 



0.15 



0.25 



0.35 



0.45 



FIG. 4: The generalized Stokes-Einstein-Debye relation between Dr{x,<j)) and high-frequency vis- 
cosity riooix,(j)) is not satisfied. Solid lines are interpolating spline fits to our simulation results 
(symbols). 
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The GSED relation for rotational diffusion is more strongly violated than its translational 
counterpart. As shown in (sil, {Dt/D^) x (?7oo/'7o) increases practically linearly in 0, even 
for non-permeable particles, whereas in Fig. H] a pronounced non-linear increase is observed. 

VI. CONCLUSIONS 

Using the HYDROMULTIPOLE simulation method, the short-time rotational self-diffusion 
coefficient, Dr{x,(f)), of uniformly permeable spheres was calculated to high precision as a 
function of permeability and volume fraction. 

An interesting scaling relation was found between Dr{x,(j)) and the corresponding coef- 
ficient, Dr{oo,(f)), of non-permeable, solid spheres of the same size, where the permeability 
enters only through the first-order rotational virial coefficient. A similar scaling was found 
for translational self-diffusion. 

The combination of the scaling relations with accurate 2nd-order and 3rd-order concen- 
tration expansion results in Eqs. (11711181) for Z}.,.(oo,0) and Dt{oo,(f)), respectively, has led 
us to the expressions 



for the permeability- dependent self- diffusion coefficients. In combination with table |TT] for 
Xr and Xt, these are convenient expressions useful in diffusion measurement analysis of 
permeable particle systems, and as input to theories of long-time dynamic properties. The 
accuracy of these expressions is better than 1.5% for rotational and 3.5% for translational 
self- diffusion, for the whole range of volume fractions (p < 0.45 provided x > 5. We expect the 
expressions to be useful in the experimental analysis of self-diffusion, to gain a quick estimate 
of the mean porosity in concentrated systems. Moreover, they can serve as short-time 
inputs into theoretical methods of calculating frequency-dependent and long-time diffusion 
properties, such as in mode-coupling and dynamic density functional theory methods. 

The simulation results for Dr{x, 0), and recent results for rj^{x, (p), were used to scrutinize 
the validity of a generalized Stokes-Einstein-Debye relation in its dependence on permeabil- 
ity. We found this relation to be significantly violated for non-dilute suspensions, unless 
the permeability is unrealistically large. The GSED test for porous particles presented in 




(21) 



(20) 
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this paper complements earlier GSE performance tests 3l|] of different short-time diffusion 
properties. Of all considered GSE relations, only the one for the cage diffusion coefficient 
can claim a certain validity when applied to neutral porous particles [31]. However, also this 
relation becomes invalid when the particles are significantly charged |7|. 

With the present paper on self-diffusion in combination with earlier simulation results for 
other dynamic properties such as the hydrodynamic function and viscosity, and the develop- 



ment of accurate analytic approximations for these properties [8|,|31 33], we have obtained an 
essentially complete description of the short-time dynamics of uniformly permeable particles 
with no-overlap interactions. 

Dispersions of spherical particles with more complex internal hydrodynamic structure, 
such as core-shell particles, and different direct interactions, will be the subject of a future 
study. 
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